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Abstract: The nucleation in the Ising ferromagnet has been studied by 
Monte Carlo simulation. Here, the magnetic field is spreading over the space 
in time. The nucleation time is observed to increase as compared to that in 
the case of static field. The cluster of negative spins is observed to grow from 
the center. The growth of effective magnetisation is studied as a function of 
the temperature and the strength of spreading magnetic field. The ratio of 
nucleation time and effective time is also studied as a function of the strength 
of spreading magnetic field. The effective time would introduce itself as a new 
scale of time in the case of nucleation by spatially spreading magnetic field. 
The effective time and the nucleation time both are observed to decrease as 
a power law fashion with the spreading rate of magnetic field. 
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I. Introduction: 

The dynamical aspect of the Ising ferromagnet, particularly concerned 
with the nucleation and growth of clusters, has become an active area of the- 
oretical research over a very long period of time[l]. The lifetime of metastable 
states and its dependences on the temperature and applied magnetic field are 
main focus of attention. Extensive Monte Carlo studies were performed in 
last few decades[2, 3, 4]. The simulational results show strong agreement 
with classical nucleation theory. Some of the important studies may be men- 
tioned briefly as follows: (a) classical nucleation theory was reexamined by 
extensive Monte Carlo simulations and connection with hysteresis was found 
[2] . In this case, the applied magnetic field was taken uniform over the space 
and constant in time, (b) Later, nucleation in the Ising ferromagnet was 
studied in presence of a time dependent magnetic field to analyse the hys- 
teresis phenomenon. However, the field was uniformly distributed over the 
space [1]. The magnetisation reversal was studied in the Ising ferromagnet 
by a pulsed field [5]. The magnetisation reversal was also studied in the Ising 
ferromagnet by a periodic impulsive field [6]. (c) Magnetisation switching 
behaviour was also studied [4] in vector spin model with long range dipo- 
lar interaction by Monte Carlo simulation. The magnetisation reversal and 
switching were studied experimentally [7, 8]. It may be noted here that all 
theoretical and experimental investigations were made in the presence of ei- 
ther a static magnetic field or by a field varying in time [3]. But in the case 
of time dependent magnetic field the field was distributed uniformly over 
the space. No such study was found where the spatio-temporal variation of 
field was considered. It is interesting to know how the magnetsation reversal 
happens when the field is spreading spatially in time. 

In the present paper, the nucleation in the Ising ferromagnet is studied 
where the field is spreading spatially in time. The Monte Carlo simulation 
with Metropolis single spin flip dynamics is used here. In section-II the 
model and the simulation scheme are described. The simulational results are 
reported in section-Ill. The paper ends with a summary in section-IV. 

II. The model and simulations: 

The Hamiltonian representing the Ising ferromagnet (described on a square 
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lattice of linear size L) in the presence of a magnetic field is 

H = -JJ2S!S*-Y,h(r,t)S?. (1) 

<ij> i 

Where, S?(= ±1) is the Ising spins, J is the ferromagnetic (J > 0) interaction 
strength and h(r, t) is the space and time dependent magnetic field. The 
spatially spreading magnetic field is taken as 

The field h(r, t) is acting along -z direction and spreading in time (t) radially 
(r 2 = (x — x ) 2 + (y — Ho) 2 )- This is like a mexican hat spreading radially 
keeping its height (ho) fixed. The factor ho is defined as the strength of the 
field. The value of the field decreases from that at the origin (x , yo) as one 
go radially away from the origin. The effective radius (V e //) is defined as 
that value of the radius where the value of the field becomes - times the 

e 

value at origin. So, r e // = (2(a 2 + b 2 t 2 ). At time t — the effective radius 

was ay/2. As time passes the effective radius r e // increases and eventually 
the field spreads over the whole lattice. The nearest neighbour ferromagnetic 
uniform (J = 1) interaction is considered here. The boundary condition is 
periodic in both sides. 

A square of linear size L = 101 is considered and the origin of spreading 
field h(r,t) is taken at the center (x = 51, y = 51) of the lattice. Initially 
(at t — 0) value of all spins (Sf) are taken up (+1). In the random spin 
updating scheme, a site is chosen randomly and the energy of the system is 
calculated (by using equation (1)). The probability of flip of this chosen spin 
is obtained from the Metropolis rate [9] 

A H 

P(S> -St) = Min[exp(-— ), 1] (3) 

where T is the temperature, k is Boltzmann constant and AH is the change 
in energy due to spin flip. Using Monte Carlo method, this spin flipped 
(Sf — > —Sf) if the value of the Metropolis rate is less than or equal to a 
fraction (0 < / < 1) chosen randomly. L 2 such random updates of spins 
constitute one Monte Carlo Step per Spin (MCSS) and act as the unit of 
time in the simulation. 
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III. Results: 

The Monte Carlo simulation was performed on a two dimensional lattice 
of linear size L. Here L = 101, so that the site (51,51) becomes the central 
site. Initially the spreading field has a width determined by the variable 
a = 0.1 and was kept constant throughout the simulational study. Later, 
the field was allowed to spread according to the equation (2) and the spins 
were updated (randomly) with probability (given in equation -3). The total 

s z 

magnetisation m(t) is calculated as m(t) = ~y and if its value becomes 
less than or equal to 90 percent of its initial value , it nucleates and the 
time required for this is called the nucleation time r nucl . This prescription of 
calculating the nucleation time was used earlier [2] 

It was observed that a single compact cluster of down spins embedded 
in up spins grows from the central site. A typical growth of a cluster of 
down spins is shown in figure 1. This resuls is quite obvious here, since the 
field is spreading radially. As time passes the boundary sites of negative spin 
cluster experiences larger values of negative field and get flipped. In this way 
the cluster of down spins grows. From figure 1(a) and figure 1(b), it is also 
observed that, in earlier time, the topology of the lattice is not reflected in 
the growing cluster of negative spins. However, for later time, as the field 
spreads radially, this topological symmetry of the lattice starts to appear 
from the shape of the growing cluster. It may be noted that Figure-l(c) and 
figure 1(d) show the cluster of nearly square shaped. 

The growth and the nucleation studied here (in the presence of a magnetic 
field spatially spreading in time) is compared with that in the presence of 
a uniform and static magnetic field. The nucleation time was calculated as 
averaged over 20 different samples. In this case one expects that for higher 
values of field strength the nucleation time (for spreading field) would be 
many times higher than that for uniform and static field. The reason is very 
clear here. In this case all spins will not experience the maximum amount of 
field at any given instant (MCSS). However, this time (for spreading field) 
would become of the same order of that for uniform and static field for the 
low fields. In this case the nucleation time is quite large for the case of static 
field. And for spreading field, the field will get sufficient time to spread 
over the whole lattice. As a result, the nucleation time for spreading field 
will be of the same order of that of static field. These results are depicted 
in figure 2. This study clearly indicates the difference in time of the growth 
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and nucleation for different types of fields (uniform-static type and spreading 
type). 

The total number of downspins, N±, within the circle of effective radius 
(r e ff) is also calculted. Hence, the effective magnetisation, m eS (= N 1 ), is 
calculated. Here also m e s is calculated as averaged over 20 different samples. 
As time passes, the m e g increases. The growth of effective magnetisation 
(m e ff) at a fixed temperature (T = 0.85) and for different values of strength 
of magnetic field (h = 0.85, 0.80 and 0.70) are studied. The values of a 
and b in the spreading field are taken here 0.1 and 1.0 respectively. These 
are shown in figure 3. For any fixed value of field strength, m e s first grows 
very slowly and then suddenly it starts to grow very rapidly. A new time 
scale, called effective time (r e fr), may be defined here. The time required by 
the effective magnetisation to reach the value 0.01 is defined as the effective 
time T e ff . Naturally, the r e g- will increase for smaller values of the strength of 
spreading magnetic field. This feature is clear from figure 3. 

The growth of m e s is also studied for a fixed value of field strength (ho = 
0.85) and different values of temperature (T=0.90, 0.80 and 0.75). The values 
of a and b in the spreading field are taken here 0.1 and 1.0 respectively. This 
is shown in figure 4. Here, also it is observed that the effective time (r cff ) 
increases as the temperature decreases. 

To consider the effecive time (r e s) as a new time scale, in the case of 
nucleation in the Ising ferro magnet, would be justified in the following study. 
This is the main part of this paper. The ratio (R) of nucleation time (r nucl ) 
to the effective time (r eff ) is studied as a function of the strength of spreading 
magnetic field for different values of temperature for fixed a = 1.0 and b = 
0.1. This is shown in figure 5. The ratio (R) is observed to be a nonmonotonic 
function of the strength of spreading magnetic field. As the field increases 
the ratio increases first, reaches a maximum, and then it starts to decrease. 
This indicates that the effective time (r eS ) may be considered as to define a 
new time scale in the problem of nucleation in the Ising ferromagnet in the 
presence of a spreading magnetic field. 

In figure 6, the R is plotted against ho for different system sizes (L =51, 
101 and 201) at fixed temperature T = 1.0. Here, a = 1.0 and b = 0.1. 
This is shown in figure 6. It is observed that as the system size increases 
the maximum of R becomes sharper and shifts towards the lower values of 
the strength of the spreading field. However, the nonmonotonic variation 
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remains unaltered leading to the existence of new time scale in the case of 
nucleation in the presence of spreading magneic field. 

The maximum value of the ratio (-R max and the strength of the field 
(/i™ ax ) for which R becomes maximum, is plotted against the temperature 
(T) and shown in figure 7. Both /i™ ax and R max decreases as the temperature 
increases. 

The effective time (r e fr) and nucleation time (r nuc i) are studied as a func- 
tion of the spreading rate (b) of magnetic field and shown in figure 8. The 
preliminary results show that both times behave like fr -0 - 55 . 

IV. Summary: 

In this paper, the nucleation in the Ising ferromagnet is studied by Monte 
Carlo simulation using the Metropolis single spin flip dynamics. All earlier 
studies were performed with magnetic field uniformly distributed over the 
space and constant in time. A few studies are done by using a slowly varying 
(in time) magnetic field but uniform over the space. However, the present 
study differs significantly with all other earlier studies. Here, the applied 
magnetic field has a spatio-temporal variation. The applied magnetic field is 
spreading over the space in time. A Gaussian-like simple form of spreading 
magnetic field is taken here. The width of the Gaussian is spreading in time 
keeping the altitute fixed. 

The nucleation time would become larger (for fixed values of temperature 
and field strength) in comparison to that observed in the case of static and 
uniform magnetic field. In the case of static and uniform magnetic field, the 
nucleation time provides the scale of time. A considerable number of studies 
are done using this time scale which was also considered as the lifetime of the 
metastable state. However, in the case of spreading magnetic field, apart from 
the nucleation time, a new scale of time appears. This is the effective time. 
If the effective time time would not provide another time scale, this would 
be a simple factor of nucleation time. And if the ratio of nucleation time 
and the effective time was plotted against the applied field strength, it would 
show a straight line parallel to the horizontal axis. However, interestingly it 
is observed in the present study the ratio as a function of the field strength, is 
nonmonotonic. Hence, one would say that the effective time and nucleation 
time provide separately different scales of time in this problem of nucleation 
in spreading magnetic field. The existence of new time scale in the case of 
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nucleation in the Ising ferromagnet is reexamined for different system sizes. 

The effective time and the nucleation time are studied here as a function 
of the factor predominantly determining the rate of spreading of the magnetic 
field. It is observed that both show a power law type decay. 

The study of the nucleation in the presence of a spreading magnetic field 
would help to sustain the metastability for a longer period if one uses the 
magnetic field spreading very slowly. In technology, if slow switching is nec- 
essary one may employ this technique. One may tune the rate of magnetic 
switching by varying the rate of spreading of magnetic field. Hopefully, this 
study will also explore some new class of problems in nonequilibrium statis- 
tical mechanics in near future. 
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Fig.l Growth of negative spin clusters (on a 101 x 101 square lattice) with 
time (t). Here, T = 0.8 and h = 1.5. (a) t = 100 MCSS, (b) t = 150 MCSS, 
(c) t = 200 MCSS and (d) t = 250 MCSS. 
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Fig. 2 The nucleation time (r ra ) is plotted against l/Zio, for uniform-static 
field (+) and spatially spreading (x) field. Here the temperature T = 0.8. 
For spreading field a = 0.1 and b = 1.0. 
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Fig. 3 The effective magnetisation (m e s) plotted against time (t). Different 
symbols corresponds to different values of field strength (ho)- ho = 0.85(+), 
h = 0.80(x) and h = 0.70(*). Here, the temperature T = 0.85 and a = 0.1 
and b = 1.0. 
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Fig. 4 The effective magnetisation (m e s) plotted against time (t). Different 
symbols corresponds to different values of temperatures (T). T = 0.90(+), 
T = 0.80(x) and T = 0.75(*). Here, the value of field strength h = 0.85 
and a = 0.1 and b = 1.0. 
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Fig. 5 The ratio (R = ^md.) j s plotted against the field strength (ho) for 
different temperatures (T). Different symbols corresponds to different tem- 
peratures. T = 1.3(D), T = 1.2(A), T = 1.0(+), T = 0.8(x), T = 0.7(o) 
and T = 0.6(*). Here, a = 1.0 and b = 0.1. 
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Fig. 6. The ratio R(— is plotted against field ho at fixed temperature 
(T = 1.0) for different system sizes. Here, a = 1.0 and b = 0.1. Different 
symbols represent different sizes, L = 51(+), L = 101(x) and L = 201(*). 



14 



3.5 



2.5 



n 1 1 1 1 1 r 



* • 

* * f 

J I I I I I L 



1.5 

0.2 0.4 0.6 0.8 1 1.2 1.4 
T 

Fig. 7 The /i™ ax («) and -R ma x(*) plotted against the temperature (T). 
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Fig. 8 The nucleation time r nuc i (x) and the effective time r e fr (+) are plotted 
against the spreading rate (b) for the magnetic field in logscales. The solid 
line is y = 50a;~ - 55 plotted in log-log scale. 
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